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Abstract 

By studying a modified Friedmann equation which arises in an extension of gen- 
eral relativity which accommodates a time-dependent fundamental length L(t), we 
consider cosmological models where the scale factor diverges with an essential singu- 
larity at a finite future time. Such models have no dark energy in the conventional 
sense of energy possessing a truly simple pressure-energy relationship. Data on su- 
pernovae restrict the time from the present until the Rip to be generically longer 
than the current age of the Universe. 
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1 Introduction 



The cosmic concordance of data from three disparate sources: Cosmic Microwave Back- 
ground (CMB), Large Scale Structure (LSS) and Type la Supernovae (SNela) suggests 
that the present values of the dark energy and matter components, in terms of the critical 
density, are approximately fix — 0.7 and Q m ~ 0.3. 

The equation of state of the dark energy w = p/p suggests the possibility that w < — 1, 
first studied by CaldwelljT] and subsequently in a number of papers[2l Bj- 

The conclusion about the make-up of our Universe depends on assuming that general 
relativity(GR) is applicable at the largest cosmological scales. Although there is good 
evidence for GR at Solar-System scales [4 there is no independent evidence for GR at 
scales comparable to the radius of the visible Universe. This involves an extrapolation in 
scale comparable to that from the weak to the grand unified scale in particle physics (some 
13 to 25 orders of magnitude) and could be called the cosmological desert hypothesis. The 
expansion rate of the Universe, including the present accelerating rate of cosmic expansion 
can be parameterized in the right-hand side of the Friedmann equation by including a 
dark energy density term with some assumed time dependence on the scale parameter 
a(t): pz>E ~ aP with (3 = —3(1 + w). This is a rather restricted function if we assume that 
the equation of state w is time-independent. But as soon as we admit that it may depend 
on time w(t) then the function on the right hand side of the Friedmann equation becomes 
completely arbitrary, just as does its designation as a dark energy density. We restrict the 
designation "dark energy" for the normal case where general relativity is assumed at all 
length scales and the equation of state w = p/p is a constant or at most a very simple 
function of time. 

Present data are fully consistent with constant w = —1 corresponding to a cosmological 
constant which we may accommodate on the left hand side of the Friedmann equation 
describing the expansion rate of the Universe. But cases with w ^ — 1, including w < — 1, 
are still permitted by observations. In this case, there is a choice between cooking up a 
"dark energy" density with a particular time dependence Pde ^ al3 ^ on the right-hand 
side of the Friedmann equation or changing the left-hand-side by changing the relationship 
between the geometry and the matter density, i.e. by changing GR. 

The case constant w < — 1 has the interesting outcome for the future of the Universe 
that it will end in a finite time at a "Big Rip" before which all structure disintegrates j^]. 

In the present article, we shall study an amalgam of the modification of GR due 
to Dvali, Gabadadze and Porrati[U] (DGP) and the idea of a Big Rip, in fact here a 
Bigger Rip. The aim is mere;y to study the range of possibilities in modified Friedmann 
equations but the results are sufficiently interesting to examine and such modifications 
may be constrained by observational data. 
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2 Set up 



The DGP gravity[6j arises from considering the four- dimensional gravity which arises from 
five-dimensional general relativity confined to a brane with three space dimensions. The 
underlying action is: 

S = M{tf J d 5 xVGK® + M 2 Planck J d'x^R (1) 

where and R are the scalar curvature in 5- and 4- dimensional spacetime respectively, 
and G and g are the determinant of the 5- and 4- dimensional spacetime metric. 

This leads to an interesting modification of GR which embodies a time-dependent 
length scale L(t) = Mp lanck /M(t) 3 . For cosmology it is natural to identify, within a 
coefficient of order one, L(t) with the Hubble length L(t) = H(i)~ l at any cosmological 
time t. Actually we are slightly generalizing the original DGP approach to include time 
dependence of the length scale L(t). 

Taking the four dimensional coordinates to be labeled by i, k = 0, 1, 2, 3 leads to the 
following modification of Einstein's equation, as a result of varying the action Eq.(jTj): 

where the specialized notation [()] in Eq.© means the following: 

j dx[(f(x))] = f'(0)5(x). (3) 

and in Eq.([3*|) we here identify x with the additional space dimension. 

It is interesting to generalize the Schwarzschild solution to this modification of GR[7]. 
One finds that the modification of the Newton potential at short distances is given by: 
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where r g = 2Gm is the Schwarzschild radius. 



The fractional change in the Newtonian gravitational potential at cosmological time t 
at orbital distance r from an object with Schwarzschild radius r g is therefore 
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L(t)\ 



(5) 
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In the Bigger Rip scenario we will describe the characteristic length L(t) will decrease 
with time according to 

L(t) = L(t )T(tr (6) 
where the power satisfies p > 1 (p < 1 implies that L(t) would increase) and where 



T(t) 



(trip t) 
(trip ^o) 



(7) 



in which t r i v is the time of the Rip. A bound system will become unbound at a time tjj 
when the correction to the Newtonian potential becomes large. We make adopt the value 
of tu defined from Eq.© by 



We can rewrite Eq.(jHJ) as: 
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where 7 = (t rip -t ) 1 . 

We shall define another later time t caus as the time after which the two objects of a 
bound system become causally disconnected from t caus until t r i p . This is defined by the 
equation: 

Iq 1 a(t caut 
c \ a(tjj) 



(t 



rip 



tr 



(10) 



As an example taking p = 1 with the values Lq = H 



(UGy) 



-1 



1.3 x l^cm 



and 7 = (2§Gy) 1 we arrive at the entries in the following Table: 



Bound system 


k(cm) 


r 9 (cm) 


(trip tu) 


(trip t caus ) 


Typical galaxy 


5 x 10 22 


3 x 10 16 


lOOMy 


4My 


Sun-Earth 


1.5 x 10 13 


2.95 x 10 5 


2mos 


31hr 


Earth-Moon 


3.5 x 10 1U 


0.866 


2weeks 


lhr 



Note that the values we find for (t r i P — tu) are consistent with those found in [Hj. The 
corresponding dark energy would have equation of state w = — 1 — I7L0 = —1.466 which, 
like that of jB], is now outside of the range allowed by observations [8 J if we assume a 
constant equation of state although it is allowed in the present model with its time- 
dependence. As another example, with more normal present w, we can increase the time 
to the Rip to 7 = (50G?/) -1 in which case w(t ) = —1.19 and the Table is modified to: 
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Bound system 


l (cm) 


r g (cm) 


(trip t{J ) 


(t r ip t caus j 


Typical galaxy 


5 x 10 22 


3 x 10 16 


250My 


7My 


Sun-Earth 


1.5 x 10 13 


2.95 x 10 5 


5mos 


2days 


Earth-Moon 


3.5 x 10 1U 


0.866 


lmo 


2hrs 



so with the more lengthy wait until the Big Rip the disintegration of structure and causal 
disconnection occur correspondingly earlier before the eventual Rip. 



3 The Bigger Rip 



The modified Friedmann equation for DGP gravity is 
so that we arrive at: 

- a =H(t) = H{U)±- p (12) 



hanck) 



In Eqs. (jllll2j) we can neglect, for the future evolution, the term (p M + p 7 )/(3ikf 
on the right- hand- side of the modified Friedmann equation. Defining 7 = —dT/dt 
(Uip - to)' 1 gives: 

.Tit) dT 



and hence, for p = 1, which is similar to dark energy with a constant w < — 1 equation of 
state: 

a(t) = T~^k) (14) 
while for the Bigger Rip case p > 1 one finds 



a(t) = a (t ) exp 



1 1 



Tf- 1 ; (p-l) 7 L(* ) 



(15) 



Here we see that the scale factor diverges more singularly in T for p > 1, hence the 
designation of Bigger Rip. In particular we study the values p = 2, 3, • • • as alternative to 
the "dark energy" case p = 1. 
Inverting Eq.(|15J) gives: 

T=[l + (p-l)jL(t )kia(t)]-tt (16) 
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In this case there is strictly no dark energy, certainly not with a constant equation of state, 
but we can mimic it with a fictitious energy density pl by noticing that H 2 ~ T~ 2p and 
writing 

PL ~ [1 + (p - lhL(t ) In a(t)] (17) 
If we use Eqs. fll5|) and (|T7jl in conservation of energy 

1(^3) = _ p | (a 3) = -^(^l^) (18) 

we find a time-dependent u>i(t) for the "fictitious" dark energy 

2 dL(t) , \ 

«,,(*) = -l + -^ii (19) 

= _1_2 ?7^o) 

31 + (p-l) 7 L(t )lna(t) V ; 

so the effective wi(t) has the limiting values Wi(to) = — 1 — |p(7L(t )) and WL(t rip ) = — 1. 

We may check consistency with the space-space components of Einstein's equations 
which are 

£ = -i*!-47rGp (21) 
a 2a z 

which leads to 

H = -47rG(p L + p L ) = -Anp L (l + w L ) (22) 
so that with H = L~ x and p^ = 3H/(87rGL) we find a w^it) consistent with Eq. (J2(J|) . 



Keeping the value L = Hq 1 = (UGy)- 1 = 1.3 x 10 28 cm and putting 7 = (20Gy)~ 1 
and p = 2 we arrive at the entries in the following Table: 



Bound system 


l (cm) 


r 9 (cm) 


(trip tjj ) 


(t C aus tjj ) 


Typical galaxy 


5 x 10 22 


3 x 10 16 


2.37Gy 


1.14Gy 


Sun-Earth 


1.5 x 10 13 


2.95 x 10 5 


9.6 x W 4 y 


7y- 


Earth-Moon 


3.5 x 10 10 


0.866 


2.5 x 10 A y 


6mos. 



Note that for the p = 2 case we have tabulated the difference (t caus — tu) rather than 
(t r i P — tcaus) because in this case the expansion is so rapid. 



Next we turn to observational constraints on the parameters Lq and 7 for p = 2. 
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4 Observational Constraints 



In the previous section, we discussed the future universe in the model. In this section, we 
discuss constraints on model parameters from SNelA data. To discuss the constraint, we 
have to include other component such as cold dark matter (CDM) and baryon. Including 
all components, we can write the Friedmann equation as jUj 



Pm 


1 


3^Planck 


f 4L 2 


= Pm/ Pcrit 


— PmO 




If we define the density parameter fl m = p m j 'p crit = p m o(l + z) 3 , we can rewrite Eq. 



as 

where Qk and Ql are defined as 



H 2 = Hi 



(24) 



■k 



Thus at the present time, we have the relation among the density parameters, 

tt k + (yvT L + y/n L + Qm) 2 = 1. (26) 

To obtain a constraint from the SNela observations, we assume that red-shift depen- 
dence of Ql is as in Eq. (JT7jl and that L is dependent on time as Eq. (|B|); thus we can say 
that we consider a new component p^ which is defined as Eq. (|17|). 

Now we discuss the constraint on this model from SNela data using recent result jH]. 
In Fig. [TJ we show contours of 95 and 99 % C.L. in flL-fl m plane for jLq = (which is 
the constant L case), 0.5 and 1 with p = 2. In the figure, we also plot the line for the 
flat universe. Notice that the line is different from the standard case because we have the 
modified Friedmann equation in this model. To obtain the constraint, we marginalize the 
Hubble parameter dependence by minimizing x 2 f° r the fit. 



In Fig. |2l we show the constraint on Q m - , yLQ plane assuming the flat universe (Qk = 0). 
If we take the value fi m = 0.3, we can find an upper limit for ^L which is 7-L ^ 0.7. This 
implies that the time remaining from now until the Rip is constrained to be generically 
at least somewhat longer than the current age. If we make L larger than the length 
corresponding to the age of the Universe then the upper bound on 7 diminishes and hence 
the time until the Rip increases. 

We note that because the effective equation of state w(t) is varying with time its present 
value w(to) can be more negative than allowed by constraints derived from assuming 
constant w. Our constraint on < 0.7 permits w(t ) = —1 — |p7-L to be as negative 
as w(to) = —1.9 for p — 2. Assuming constant w, on the other hand, gives jH] w > —1.2. 



6 



5 Discussion 



The present article is a natural sequel to our previous paper jS] about dark energy in which 
it was pointed out that no amount of observational data can, by itself, tell us the fate of 
dark energy if we allow for an arbitrarily varying equation of state. The three possibilities 
listed were: there may firstly be a Big Rip, or secondly dark energy may dominate but with 
an infinite lifetime or thirdly the dark energy may eventually disappear leaving a matter- 
dominated Universe. Given that observational data are insufficient, only a successful and 
convincing theory may inform us confidently of the future of the Universe. 

The Big Rip was the most exotic of the fates and there seemed tied to a phantom 
w < — 1 dark energy. However, here we have studied a Bigger Rip, in which the scale 
factor is even more divergent at a future finite time than for the Big Rip, which is achieved 
by modifying gravity and omitting dark energy. In the model, as with the phantom 
case, structures become unbound and subsequently their components become causally 
disconnected before the Universe is torn apart in the Rip. 

Work by other groups ^U] has recently suggested that quite different values of cos- 
mic parameters can be acceptable if one relaxes the most conventional and conservative 
assumptions of general relativity at all length scales and a dark energy with a constant 
equation of state. 
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Figure 1: Constraint from SNela observation in Q L -Q m plane for 7L0 = (bottom), 
7L0 = 0.5 (middle) and 7L0 = 2 (top). Contours are for 95 % (dotted line) and 99 % 
(dashed line) C.L. constraints respectively. The solid line indicates parameters which give 
a flat universe. 
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Figure 2: Constraint from SNela observation in f2 rc -f2 m plane. Contours are for 95 % 
(dotted line) and 99 % (dashed line) C.L. constraints respectively. In this figure, we assume 
a flat universe. 
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